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ABSTRACT 

Context. Abundance anomalies have been observed in field sdB stars and in nearly all Horizontal Branch (HB) stars of globular 
clusters with Tejf > 1 1 000 K whatever be the cluster metallicity. 

Aims. The aim is to determine the abundance variations to be expected in sdB stars and in HB stars of metallicities Z > 10 and what 
observed abundances teach us about hydrodynamical processes competing with atomic diffusion. 

Methods. Complete stellar evolution models, including the effects of atomic diffusion and radiative acceleration, have been computed 
from the zero age main-sequence for metallicities of Zq = 0.0001, 0.001, 0.004 and 0.02. On the HB the masses were selected to 
cover the Tetf interval from 7000 to 37000 K. Some 60 evolutionary HB models were calculated. The calculations of surface abundance 
anomalies during the horizontal branch depend on one parameter, the surface mixed mass. 

Results. For sdB stars with Teff < 37000 K and for HB stars with TefF > 1 1 000 K in all observed clusters, independent of metallicity, 
it was found that most observed abundance anomalies (even up to ~ x 200) were compatible, within error bars, with expected 
abundances. A mixed mass of ~ 10"^ Mq was determined by comparison with observations. 

Conclusions. Observations of globular cluster HB stars with Tgff > 1 1 000 K and of sdB stars with Tf-g < 37 000 K suggest that most 
observed abundance anomalies can be explained by element separation driven by radiative acceleration occuring at a mass fraction of 
~ 10"^ Mq. Mass loss or turbulence appear to limit the separation between 10"' Mq and the surface. 

Key words. Horizontal Branch ~ sdB stars - Stellar evolution - chemical composition 



1 . Astrophysical context 

Large abundance anomalies have been observed on the horizon- 
tal br anch (HB) of NGC6752 , NGC1904 NGC2808, M15 and 
M13 ('Behretal.' 1999, '2000'; Behr '2003^. 'Moehler et al.ll2000t 
iFabbi an et al. 2005; Pace et al. 2006): whereas those stars cooler 
than about llOOOK have the same composition as giants, those 
hotter than llOOOK usually have larger abundances of some 
metals by large factors. This occurs in all clusters having blue 
enough HB stars irrespective of their metallicity as defined by 
their giant branch stars. 

Field sdOB stars are observed to have lar ge abundanc e 
anomalies compared to Pop I stars ( for a review see Heber 2009). 
They were already recognized by ISargen t & Searle (1968) to 
have a surface composition different from the one with which 
they formed. While cool (Tgff < 10"* K) field Pop II stars have 
low Z, the sdB^ which must come from the red giants often 
have iron peak abundances which are solar or even larger While 
sdBs correspond to the blue end of the HB, the hotter sdO stars 
are apparently a mixed bag of post HB stars and other highly 
evolved evolutionary stages. Spectroscopically, sdOs have much 
more diverse characteristics than sdBs. In this paper we extend 
the Teff coverage of evolutionary models to sdBs but not to sdOs. 



In preceding papers (iMichaud et al.l 120071: iMichaud et al.l 

2008D the evolution of a Pop II staiQ with Zq = lO""* was fol- 
lowed from the zero age main-sequence to the middle of the HB 
and comparisons were made to observed abundance anomalies 
in M 15. The overabundances are explained by atomic diffusion 
driven by radiative accelerations in stars with > 1 lOOOK and 
the sudden break in anomalies at 1 lO OOK was shown t o be re- 
lated to observed rotation velocities (Ouie vv et al.ll20 09\. Given 
the relatively large observational error bars, the anomalies ap- 
peared compatible with a simple diffusion model involving only 
one parameter, the m ass of the out er region mixed by turbulence. 
It was determined bv lMichaud et al. (2008) to be about 10"^ M,. 

Extending these calculations to higher metallicities allows 
comparing to field sdOB stars. While these have the disadvan- 
tage that their original metallicity is unknown, some of them are 
much closer than any globular cluster and so their surface com- 
position can be determined more precisely than that of corre- 
sponding HB stars of clusters. This opens the possibility of fur- 
ther constraining the process competing with atomic diffusion in 
causing abundance anomalies. 

According to Mo rrison et al.l (|2009i) and iKinman et al.l 
( 2009) . inner halo stars have a mean [Fe/H] = -1.6 and almost 
all have [Fe/H] < -0.8. Presumably most field sdB stars have a 



' We use the expression HB stars for Horizontal Branch stars in clus- 
ters, and sdB stars for those in the field. In practice, since calculations 
are always done for a given metallicity, they are always done for cluster 
stars causing apparent contradictions to that rule when comparisons are 
made to field stars. 



- In this series of papers, the relative values of the a elements are 
increased following VandenBerg et al. (2000). Models are labeled ac- 
cording to their or iginal Zq value calc ulated before the a coiTection. 
See also Table 1 of lRichard et al.l2002l . This correction was not applied 
in the Zq = 0.02 models. 
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similar range of original metallicities. It is then useful to deter- 
mine to what extent the surface abundances of sdBs should be 
affected by the various metallicities they may have formed with. 
Carrying calculations with different metallicities allows not only 
comparisons to globular clusters whose giant branch shows vari- 
ous metallicities but also to analyze potential causes of the abun- 
dance range observed in field sdBs. 

Among previous work relevant to this paper, one may men- 
tion that radi ative accelerations of m etal s and He in sdBs wer e 
calculated by iBergeron etaH (Il988b and iMichaud et aH (Il989h . 
The calculations of ^rad(He) were for stars of 40000 K or more. 
The role o f He diffu s ion fo r the structure of sdBs has been inves- 
tigated bv lHu et al.l(l2009h and lHuetall f2010) but without the 
important contribution of metal diffusion with radiative acceler- 
ations. 

The sdBs correspond to the highest Teff HB stars observed 
in globular clusters. Why did they lose more mass above the He 
burning core than the cooler HB stars? A number of more or less 
complicated scenarios have been suggested for them. Here we 
assume a simple scenario and determine what anomalies would, 
in this scenario, be expected on the surface. We do not claim to 
exclude more complicated ones but consider useful to determine 
if surface abundances can be understood in a simple scenario, 
assuming nature chose that way. 

In this paper, the evolution is carried out with models evolved 
with atomic diffusion from the zero age main -sequence and 
thro ugh the Red Giant Branch, as described in iMichaud et al.l 
(1201 0). During the HB, it is continued with a surface mixed zone 
with enough turbulence to force abundance homogeneity. The 
mass of this mixed zone is kept constant. Abundance anoma- 
lies depend on its extent which is determined by a comparison 
to observations (see Sect. l3.2l i. In general the calculated anoma- 
lies also depend on the original metallicity of the star. In the 
case of globular clusters the latter is fixed by the metallicity of 
the giants of the cluster Th is approach is different from that of 
ICharpinet et all ( [1997) and iFontaine et"ari ( .2003 ), in that these 
authors assumed the abundances in the exterior regions to be 
such that gravity and radiative accelerations were at equilibrium. 
This led to a parameter free determination of abundances which 
is a function of T^ff only. Their model has considerable success 
in explaining observations and pulsation properties. However 
can the assumed equilibrium always be reached during evolu- 
tion on the HB? Is there a sufficient Fe reservoir to fill the re- 
gion of the Fe opacity bump? Does a mixing process modify it? 
Observations suggest a range of abundances for some species in 
sdBs at a given Teff. 

After a very brief description of the calculations (Sect.|2), 
the internal structure is analyzed (Sect.|3), contrasting that of 
14000 and 30100 K stars; the resulting surface abundances are 
discussed and compared to observations in Sect.|4] insisting on 
the metallicity dependence of clusters which leads us to continue 
with field sdB stars (Sect. l4.4] i. After a summary of the main re- 
sults (Sect. l5.ll l, the potential role of mass loss (Sect. l5.2T i and 
jU gradient inversion (Sect. 15.3^ are mentioned. While this paper 
is concerned with the role of stellar evolution for surface abun- 
dances, it is concluded by a brief discussion of the need to test 
the models using asterosismology (Sect. l5.4] i. 

2. Calculations 

Stellar evolution models were calculated from the zero age 
main-sequence to the HB as described in Michau d et al. (2007). 
Using opacity spectra from llglesias & RogersI (IT996I) all as- 
pects of atomic diffusion transport are treated in detail from 
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Fig. 1 Turbulent diffusion coefficient in a 0.50 Mq model. The 
left axis is for the diffusion coefficient and the right axis for 
the Fe mass fraction. The coefficient is anchored at a mass 
AMo = 10"^^ Mq. The vertical dashed line is at AMq. The Fe 
mass fraction, X(Fe), is nearly constant from the surface down 
to AM - 10"^ Mq which corresponds to Fe being mixed down 
to approximately 3 x AMq. 

fir st principles. These m odels are called models with diffusion 
in lMichaud et alj(l2007l) . 

During the HB, the surface convection zone includes very 
little mass for T^ff > 10000 K. To carry out HB evolution, it 
was found necessary to keep a small surface zone mixed. A sim- 
ple turbulent diffusion coefficient (see Fig.[TJ, similar to those 
used for AmFm stars of the main-sequence, was used to mix 
the exterior region during HB evolution. Turbulence has been 
assumed large enough to mix completely the regions between 
superficial convection zones; suc h mixing is expected from the 
results of numerical sim ulations (iKupka & Montgomervll2002l; 
'Frevtag & Steffen"2004'). Below the deepest surface convection 
zone, the turbulent diffusion coefficient has been assumed to 
obey a simple algebraic dependence on density given, in most 
calculations presented in this paper, by 



Dt = 10'*D(He)o'^° 



(1) 



where D(He)o is the atomic diffusion coe fficient of Hqj at som e 
reference depth. For most calculations in lMichaud et al.l (l2008h . 
the turbulent diffusion coefficient was anchored at a given con- 
stant temperature, Tq. Then 



Po = piTo) 



(2) 



and equation (|2]l is given by the stellar model. That density varies 
during evolution. While a given po or To during the evolution of 
one star of a given mass on the HB, approximately occurs at a 
constant value of AM = M» -M,-, the latter changes considerably 



' The values of £)(He)o actually used in this for- 
mula were always obtained — for programming con- 
venience — from the simple analytical approximation 
0(He) = 3.3 X lO-'^'r^ V[4pln(l + 1.125 x IQ-'^T^/p)] (in cgs 
units) for He in trace amount in an ionized hydrogen plasma. These can 
differ significantly from the more accurate values used elsewhere in the 
code. 
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Fig. 3 Mass within and above the He II, He I and H convection 
zones as a function of Tefj. The H convection zone is separate 
from the He I convection zone in only a few models. For stars 
hotter than 22000 K, the temperature at the bottom of the convec- 
tion zone is less than a factor of 2 larger than Tetf and is probably 
not determined precisely by interior models; it is not shown. The 
color code is the same as in Fig.|2] 



between cool and hot HB and especially hot sdB stars. For this 
paper, most calculations were carried out with 

po = p(AMo). (3) 

and equation (O is given by the current stellar model. In words, 
in the calculations reported in this paper, po of equation ([T]l is 
the density at which AM - AMq in the evolutionary model. For 
AMq = lO"^'^ Mq the concentration of most species is constant 
for AM ^ 10"^ " Mq when the assumed turbulent coefficient is 
10"* X the He diffusion coefficient at AMq and varying as p^*. 
As one increases AMq, one defines a one parameter family of 
models. As may be seen from Fig.lT] turbulence decreases very 
rapidly inward and using such turbulence is equivalent to mix- 
ing a certain masfl The exact mass over which concentration 
is kept uniform by turbulence depends on the driving terms in 
the atomic diffusion equation and so varies slightly with atomic 
species. As may be seen from Fig.lT] it is 3 x AMq for Fe, but it is 
10 X AMq for He (not shown). For log( AMq/ Mq) = -7.5 these 
models are labeled dM-7.5D10K-4. 



3. Abundance anomalies and internal concentration 
variations 

Calculations were carried out for metallicities of Zq =0.02, 
0.004, 0.001 and 0.0001. CoiTesponding tracks in the HR dia- 
grams and in the log g-T^ff plane are shown in Fig. |2] The time 
evolution of L, Teff, of surface abundances of He, Ca, Cr and 
Fe, and of other variables of interest may be found as an on line 
Appendix (App.lB} for all models. 

^ It was found by [Richer et alj \2OO0) that using 2 or 3 instead of 
4 as the exponent for p in equation Q only chang ed the mixed mass 
corresponding to AMq. It was further found by Tal on et al.l ( 120061) that 
distinguishing between different turbulence models required an accu- 
racy of 0.03 dex which is beyond the accuracy of current abundance 
determinations (see their Sect. 5). 



On the right panel of Fig.|2]the data from Fig. 2 of lGeier et al.l 
(1201 0*) is compared with the g vs Teff covered by our modelfl 
These cover the range of most stars in that paper for Teff < 
36 000 K. There are some lower gravity stars that have gravities 
more than a factor of 2 smaller than those of our models; they 
ar e identified usiiig a dif ferent color (red). According to Fig. 1 
of lCharpinet et al.l (12000 ;). they should coiTesppond to the end of 
the HB phase which is not treated in our mode l s. Wh en we com- 
pa re the result s to Figs. 3 and 4 of lGeier et al.l (l2010b , and Fig. 1 
of 'Geier et al.l (1200 8) we will usually identify them separately. 
Evolution was stopped when the He mass fraction dropped be- 
low ~ 60 % in the star center (after around 32 Myr on the HB) 
because the algorithm used to calculate diffusion velocities re- 
quires the presence of either He or H. For technical reasons elim- 
inating this limitation requires the rewriting of major sections of 
the code and this is outside the scope of this paper For a few 
models the calculations were stopped after less than 32 Myr (see 
App.lBll because of convergence problems. 

The mass in the surface convection zone is usually one of 
the most important stellar model properties for atomic diffusion 
since it is thoroughly mixed and separation generally occurs be- 
low. The bottom of the He II, He I and H convection zones are 
shown in Fig.|3] The bottom of the He II convection zone does 
not depend significantly on metallicity in so far as the points de- 
fine a single line (for instance, at low Teff, the magenta stars are 
on the same line as the blue stars). Jus t as occurs on the m ain- 
sequence of Pop II stars (see Fig 1 of iRichard et al.ll2005h . the 
mass in and above the surface convection zone depends nearly 
only on Teff which effectively defines the depth of the surface 
convection zone. Calculated models cooler than 22000 K are 
shown with metallicities ranging from Zo = 0.0001 to 0.02. For 
stars hotter than 22000 K, the temperature at the bottom of the 
convection zone is less than a factor of 2 larger than Teff and 
such outer structure is probably not determined precisely by in- 
terior models. It was decided not to show the mass above that 
Teff. The mass mixed by turbulence, ~ 10"^ in this paper 
(see Sect. l3.2l l, is always larger than the mass in the surface con- 
vection zones. The latter then does not play a major role. 

3.1. A 301 00 K SdB vs a UOOOKHBstar 

On Fig.m are compared the internal concentrations in a 30 100 K 
HB (0.51 Mq; right panel) to those in a 14 000 K HB star 
(0.59 Mq; left panel). Both started with the same metallicity 
(Zo = 0.0001) and have spent the same time on the HB, 25 Myr 
Given that they both started with Zq = 0.0001, the two have sim- 
ilar Fe peak surface overabundances (a factor of ~ 200 overabun- 
dance from the original abundance of the star) but very different 
He, N, Ne, Na, Mg, Al, and Si surface abundances (compare the 
black curves in the right and left panels of Fig.|5]l. However the 
lower temperature model has a radius 5.4 times larger than that 
of the hotter model. The hotter model is shown on the radius 
scale of the cooler one in the small central inset. The central re- 
gion (defined as the region where '^C has been produced by He 
burning) of the two models is quite similar By comparing the 
central region of the inset to the central region of the left panel, 
one notices that the various concentrations are about the same at 
approximately the same radius and inner mass (that is as a func- 



' We use color codes that depend on the envelope mass, that is on 
log,Q(M, - MHecore)/^o- This expression is used to define an index 
referring to 16 predefined colors. The same code is used for all figures 
involving stars with T^g > 20000 K. Another code is used when only 
stars with Teff < 20000 K are present. 




Fig. 2 HR diagram (left) for all calculated models with a mixed mass ~ 10"^ M©. Each segment represents the evolution of one 
model whose mass is color coded using a simple formula. The metallicity (Zq) of each series of models is indicated at the left of 
the sequence. The dotted part represents the first lOMyr of HB evolution while the solid line represents the rest of the calculated 
HB evolution, ususally an additional 27Myr. I n the right panel, g vs T^ff diagram of a subset of the calculations (those with T^ff > 
20000 K) corresponding to the observations of lGeier et al.ll2010l This panel contains data (gray and red dots as defined in the text) 
from their Fig. 2. 




Fig. 4 Color intensity coded concentrations in two HB stars of the same metallicity after 25Myr on the HB. Left panel with a 
Teff of 14 000 K (0.59 Mq) and right panel of 30 100 K (0.51 Mq). The radial coordinate is the radius and its scale is Hnear but 
the logarithmic value of the mass coordinate above a number of points, log AM/M,, is shown on the left of the horizontal black 
line. The concentration scale is given in the right insert. Small circles near the top of the left panel mark the extent of the surface 
convection zone while similar circles near the center of both models mark the central convection zone. The small inset in between 
the two panels shows the high Teff star, that is the right panel, on the radius scale of the low T^s star, that is the left panel. For 
-7 < log AM/Mt < -4 the concentration is quite different for many species. It is surprisingly so for C and O for log AM/M* > -2. 
See the text. A black and white version of this figure may be found in the On line Fig. lA.ll 



tion of m,, the mass within radius r) in both models. However 
while abundance variations caused by nuclear reactions extend 
over the inner 50 % of the radius in the hotter model, they ex- 
tend only over the inner 10 % of the radius in the cooler model. 
Furthermore the central region corresponds to rather different 
stellar mass fractions of the two models; it is the inner 99 % by 



mass of the 0.5 1 Mq star but it corresponds to the inner 85 % by 
mass of the 0.59 Mq model. The difference originates from the 
different peeling applied to the red giant in order to get to the dif- 
ferent Tgs of the two models: 14 % more mass has been removed 
from the red giant to get to the 30 100 K model than to the cooler 
model. In fact approximately the same inner mass is affected by 
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Fig. 5 Comparison of surface abundances expected in 14 000 (black line, left panel) and 30 100 K (black line, right panel) HB stars 
with Zq - 0.0001 whose interior is shown in Fig.|4] On the left panel, in order to compare with the turbulence model used by 
iMic haud et al. (2008) (see the text for details), the red curve is taken from the left panel of their Fig. 1 1 and the data is for B203 
from .Behr (.2003.) . On the ri ght panel the other curves correspond to stars of original metallicities of Zq - 0.001, 0.004 and 0.02 and 
the data for Feige48 is from lO' Toole & HebeJ(l2006l) . Grev is used for data points of atomic species which are not included in our 
calculations because they are not in OPAL opacities. 




log (AM/M^) 

Fig. 6 Internal distribution of Fe and Ni together with their relation to ^rad(Fe) and ^iad(Ni) in a 0.5 Mq model with Zq - 0.001 for 
four turbulence strengths, each anchored at a fixed AMq (see Eq. IT] and 131), identified on the upper panel. The solutions labeled 
-7.5, -8 and -9 are after 32 Myr on the HB while the solution labeled -10 is after 7Myr The lower panel gives /fR as a function 
of log AM/Mq. a local maximum of grad(Fe) and of kr occurs for log AM/Mq between -10 and -9 and corresponds to log T ~ 5.5. 
When AMq < 5 x IO^'^Mq the iron abundance becomes rapidly much larger than solar in contradiction with observations of sdB 
stars. The function 1/jUo is shown (solid lines) for each model in the panel above that for /c-r. In the same panel, the function 1 //ioeff 
is shown by dotted lines. It is defined in Sect. 15. 31 
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nuclear reactions in both models and the structure of the cen- 
tral core is similar in both. However the envelope has quite a 
different structure, since it extends much further in the higher 
mass model. Abundance variations caused by diffusion are seen 
to extend over the outer 40 % by radius for many species in both 
models so that, as a function of the fractional radius as used in 
Fig.|4] the outer region (log AM/M, < -5) has many similar- 
ities in the two. However between logAM/M, = -5 and the 
central region, the concentration variations are quite different. 
Furthermore even when the fractional radius is about the same 
in the two models, the thickness in physical distances differs by 
the scale factor of 5.4. 

On Fig.|5]the black curves represent the surface abundances 
of the two stars of Fig.|4l left panel, 0.59 Mq model; right panel, 
0.51 Mq model. The right panel also contains results for models 
of larger metallicities than used in Fig.|4l the original metallici- 
ties range up to Zq = 0.02. The surprise is the similarity of the 
surface Fe abundance: it is nearly solar in both the left and right 
panels and for all curves of the right panel; this corresponds to 
a factor of 1 to 200 overabundance from the original abundance 
of the star depending on the original metallicity. The Fe abun- 
dance seems to be determined by saturation of grad(Fe) where 
the separation occurs. One may relate them to the giad(Fe) at the 
bottom of the mixed zone for the turbulence model used. The 
black curves in the left and right panels are quite similar from 
Ca to Ni but differ strongly from C to K. This figure will be 
further discussed in Sect.|4] 



3.2. Determining the mixed mass 

On Fig. |6] are shown the internal distributions of Fe and Ni 
in four models of the same mass and original metallicity but 
with different turbulence. The original Fe mass fraction at the 
beginning of HB evolution was ~ 10 ' throughout the star 
with only small vari ations left over by red giant evolution (see 
lMichaudetaril2010h . The §rad(Fe) has clearly defined peaks at 
around logAM/Mg = -7.5 and at around logAM/Mo - -5 
(these may also be seen in Fig. 3 of Michaud et al. 2008). In 
between these peaks grad(Fe) remains slightly larger than grav- 
ity allowing a limited Fe flux to be pushed toward the surface 
and creating a small underabundance over the interval -5 < 
log AM/Mq < -4, the iron from that mass interval being pushed 
toward the surface. Between log AM/Mq = -5 and slightly be- 
low the bottom of the mixed zone, the Fe abundance is deter- 
mined by the requirement of flux conservation. Immediately be- 
low the mixed zone, an Fe gradient develops to satisfy the dif- 
fusion equation as §iad(Fe), the Fe abundance in the mixed zone 
and the Fe flux coming from below evolve. The largest turbu- 
lence shown is labeled by logAMo/M© - -7.5; it is used for 
most calculations and leads to about solar surface Fe abundance. 
The surface Fe abundance is mainly determined by the value of 
^rad(Fe) immediately below the mass mixed by turbulence. Since 
at log T ~ 5.4, giad(Fe) is larger than g even at solar Fe concen- 
trations, mixing only from the surface to log AM/Mq = -8.6 (or 
log r ~ 5.4) leads to larger X{Fe) than observed. One needs to 
mix further inwards, to log AM/Mq ~ -7 for surface Fe abun- 
dances to be about solar. This in practice is what fixes the mixed 
zone for the models presented in this paper The mass in the 
mixed zone was so fixed by the Fe abundance in sdBs. It was 
then used for all sdB and HB calculations with Tgff from 7000 
to 35000 K. A mixed mass a factor of 10 smaller would clearly 
lead to unacceptably large Fe abundances, while a factor of three 
larger is marginally excluded. 



In the third panel, is shown where po is the re- 

duced mass per nucleus, that is excluding electrons. When 
log(AM()/M0) < -7.5, it has a local maximum interior to the Fe 
accumulation. It is caused by the accumulation of metals pushed 
upwards by g^^. The local maximum is very small in the case 
log(AM()/M0) = -7.5. There still is an accumulation of metals 
in this case above log AM/Mq - -6.0 but its effect in increasing 
Pi) is largely canceled, for log(AMo/Mo) - -7.5, by the de- 
crease of the He concentration. It was verified (not shown) that 
in the case of turbulence with log(AMo/ Mq) - -7.5, there was 
no po inversion related to metal overabundances for metallicities 
of Zq = 0.0001, 0.001 and 0.02 and, in each case, in both a star 
of 30000K and of 15000K. The case shown in FigHis the one 
with the largest effect among all those verified. 

3.3. Radiative accelerations and interior mass fractions 

Radiative accelerations are shown at ~ 25 Myr after ZAHB in 
Fig. El in the left panel for various metallicities while in the 
right panel, they are shown at various TetfS for a metallicity of 
Z() = 0.001. The gradS clearly vary more with Tgff at a given 
metallicity than with metallicity at a given T^ff. The change with 
metallicity is mainly due to saturation of the lines which causes 
a general reduction of gi-ad as metallicity is increased; also as 
metallicity is increased, T increases at a given log(AM() / Mq) 
and a given ionization state shifts towards the surface. As T^ff 
increases however, the radiative flux increases with Tgs'^ and 
furthermore, gravity increases by orders of magnitude, strongly 
affecting the ionization equilibrium at a given log(AM()/M0) 
within the star. As seen on the figure, a given g^d maximum 
migrates toward the surface as T^ff increases. 

In the right panel of Fig.|7] one notices large scatter for 
^rad(P) in the Tgff - 10700K model. We have verified by 
an analysis of the spectra used to calculate grad(P) that very 
few lines contribute to it in the T^ff = 10700K model. At 
log(AM()/M0) - -6, where the gmdiP) curve is most irregu- 
lar, a single line contributes at least 30% of the value and 5 
lines, the following 30 %. Given that these are calculated in 
opacity sampling, this leads to the observe d fluctuations. For 
more d etails, see the disc ussion of g r ad(Li) in lRicher & MichaudI 
('2005') and in Sect. 3 of lVick et al] (|2010|) . Because of the un- 
certainties in giad(P), the corresponding curve for X(P) is not 
shown in Fig. [8] The gad& only matter where turbulence does 
not force homogenization so only for log(AMo/M0) > -7. No 
large glitches due to improper sampling occurs for any of the 
other cases shown for log(AMo/MQ) > -7. Abundance varia- 
tions as seen in Fig. [8]do lead however to variations in g,adS. In 
the same model, for instance, the Ca abundance minimum below 
the mixed zone is partly responsible for the grad(Ca) peak at the 
same place. The large overabundance of Ni in the mixed region 
of the Tgff = 10700 and 15 000 K models are largely responsible 
for the rapid variation of giad(Ni) there. 

In the left panel of Fig. [8] one sees that the surface abun- 
dances at 25 Myr increase as the starting abundances increase. 
However, saturation plays a role for all species shown: the range 
in surface abundances is smaller than the range of starting abun- 
dances. Saturation is most evident for Si and Fe, whose surface 
abundances differ by at most a factor of 3 while starting abun- 
dances covered a range of 200, still visible in the deeper layers. 
In the right panel of Fig.|8] one sees a much less regular behavior 
The starting abundances were the same for all models as still evi- 
dent from the deeper layers, at log(AM()/ Mq) = -4. The surface 
abundance of Si is smallest in the coolest model, but the sur- 
face abundance of Ni is largest in the coolest model. Iron peak 
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Fig. 7 Radiative accelerations for a subset of the included species in models with turbulence anchored at a fixed AMo - 10 '' Mq, 
after approximately 25 Myr on the HB. Dotted lines represent gravity. Left panel: in models with red ~ 3 1000 K and metallicities of 
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masses corresponding, from bottom to top, to T^^s of 10700, 15 000, 20000, 25 000 and 30 400 K. 



elements have a similar behavior with Teff. Lower mass nuclei 
vary much more differently (compare Al and Ca for instance) as 
could be expected from the behavior of g,adS on the right panel 
of Fig.|7] The shape of the surface abundance curve as a function 
of Z, the nuclear charge number, is then expected to vary with 

4. Comparison to Observations 

4.1. Metallicity, Ten, and age dependence 

Before comparing to observed abundance patterns of individual 
stars, it is useful to evaluate the dependence of calculated abun- 
dances on original metallicity, T^s and age. 

If one observes a star of 30000 K, how sensitive are its sur- 
face abundances to the original stellar metallicity? In the right 
panel of Fig. |5] the surface abundances at 25 Myr after ZAHB 
are shown in four stars of original metallicity ranging from 
Zo = 0.0001 to 0.02. They were chosen to have Teff as close 
as possible to 30000 K since Teff is usually the better determined 
quantity. This imposed different masses which range from 0.47 
to 0.51 Mq. While the original abundances vary by a factor of 
200, the final ones vary by at most a factor of 30 and for Fe and 
Ni by a factor of about 3. The largest differences are for species 
around Ca and Ne. Generally a star that starts with a larger metal- 
licity ends up having larger abundances of metals, but there are 
a few exceptions such as P and S in the Zq = 0.004 and 0.02 
models. This must be caused by a competition f or photons with, 
for instance, Fe (see Fig. 4 of lRicher et al.ll998l) . Helium is most 
underabundant in the higher Zq stars. 

Comparing the right and left panels of Fig. |5] shows that a 
factor of 2 difference in T^^ for stars originally with Zq - 0.0001 



does not change considerably the expected Fe peak abundances 
but leads to large differences for atomic elements with masses 
equal to or smaller than that of Ca. While a peak occurs for N 
in the higher T^-s model, it occurs for P in the lower T^s one. 
The differences can reach 1-2 orders of magnitude. However 
for smaller T^s differences, the abundance variations tend to be 
relatively small (see Fig.|9]i. Phosphorus is most sensitive, its 
abundance variation reaching a factor of 5 between 11500 and 
14000 K (however see also Sect. l3.3l ). But for Fe and Ni an un- 
certainty of 2500 K hardly affects the expected surface abun- 
dances even though they are more than a factor of 100 over- 
abundant. By comparing the left and right panels of Fig.fTOlone 
notes that a relatively small Teff change (~ 3000 K) can lead, in 
a cluster with Zq = 0.001, to significant abundance variations for 
species between Si and Ti. For heavier iron peak elements, the 
variations are much smaller 

Examples of age dependence are shown on Fig. [10] The 
abundance variations are rapid for the first 5 Myr but much 
smaller between 10 and 30 Myr for all species whose abundance 
increases. The variations are much more regular for species with 
A < 10 which become underabundant. In a population of stars of 
various unknown HB ages, one consequently expects most stars 
to have at least the overabundance attained after the first 5 Myr. 
See also Figs.lBT1tolB4l 

4.2. Low metallicity clusters 

In Michaud et al ] (12008 ') one turbulence model anchored at a 
fixed T was shown to reproduce abundance anomalies observed 
on the HB of M 15 from 11000 to 15000 K. A single turbulence 
model was used for all HB stars of that cluster For this paper, we 
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Fig. 9 Surface abundances expected for models with turbulence anchored at AMj Mq - 10"^^, Zq ~ 0.0001 and three slightly 
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tried using the same model for HB stars of higher metallicities 
and for sdBs of higher T^g but it was found not to reproduce all 
observations. We then searched for a model anchored at a fixed 
AM/ Mq, so at a fixed fraction of M©. It is a priori not clear how 
the different turbulence models (given by the difference between 
Eq. 12] and [3] during evolution) should affect surface abundances. 
It was first verified that one obtains with turbulence anchored 
at AM/ Mq - 10"^'^, as good an agreement bet ween observa- 
tions and calculations as shown on Fig. 9 to 12 of iMichaud et al.l 



(l2008l) bv comparing the detailed abundances for 2 stars of Ml 5 
and by comparing the Fe abundance for all stars as a function of 

Results obtained with both turbulence models, are compared 
on the left panel of Fi^.|5]to observations of the Teff ~ 14000 K 
star B203, from Behr (2003). For all calculated species, the sur- 
face abundances calculated with both turbulences are within a 
factor of 3 of each other except for S where the agreement with 
observations is better with the model of this paper Note that for 
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Fig. 1 1 Left panel Surface concentration of Fe as a function of Teff in clusters (M15, M68 and M92) from lBeh3l2003l with original 
Fe concentration ~-2.3 . Models are shown for this metallicity. Each color coded segment represents the evolution of the surface 
Fe concentration. The line is solid for the time interval from 10 to 30 Myr after ZAHB but dotted from to 10 Myr Bottom part: 
rotation velocit ies; three st ars have their rotation velocities and [Fe/H] in red to facilitate their identification. Right panel Field 
stars also from lBehrll2003[ The distribution of metallicity for the cooler field stars gives an indication of the original metallicity 
distribution for the hotter field stars as well as for sdBs. 



iron peak elements the abundances are more than lOOx the orig- 
inal abundances. 

In Fig.|9l observations (lBehdl2003h of star B315 of M15 
(Teff ~ 13 000K) are compared to calculations with turbulence 
anchored at AM/ Mq - 10"^'^ for stars of slightly different 



masses or T^ff. Two observed values are shown for Mg, the up- 
per one for Mg I and the lower one for Mg II. The agreement 
may be compa r ed t o that obtained in the right panel of Fig. 1 1 of 
lMichaudetal.l(l20"08i) . At least from their effect on surface com- 
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position, the two turbulence models are seen to be equivalent for 
B315. 

A third compari son with previous resu lts may be made by 
comparing Fig. 9 of iMichaud et al.l ( l2008l) to the left panel of 
Fig. [TT] where observed Fe abundances in the clusters M 15, 
M 68 and M 92 are compared with our calculations over the Teff 
interval from 7000 to 15 000 K. Lines represent the Fe surface 
abundance evolution during the first 30 Myr of HB evolution. 
They are color coded as a function of the mass of the model, 
the first 10 Myr, dotted, the following 20 Myr, solid. Given the 
slope of the solid lines, the Fe abundance may continue to in- 
crease by a factor of 1 .5-2 during the following 60 Myr The cal- 
culated abundances are then in agreement with observed abun- 
dances from 11 000 to 15,000K. The stars with T.ff < llO OOK 
are influenced by rotation as shown bv'Oui evv et al.l (l2009b . The 
comparison to observations is extended to more clus ters of low 
metallicity than considered by Mich aud et al.l (l2008l) . Few stars 
hotter than 1 lOOOK are added but there are three rapidly rotating 
ones; they follow the pattern explained by lOuievv et al.l (l2009t) 
and are shown in red. When compared to the results of Fig. 9 of 
[Michaud et al. (2008) those of Fig. [TT] vary less as Teff increases 
from 11000 to 16000K. There is a hint that this agrees better 
with observations, though the error bars are too large to allow a 
firm claim. 

The model used here for all calculations reproduces obser- 
vations as well as that used in Michaud et al.. (.2008) . As will be 
verified in Sect. l4.4l fitting observed anomalies of stars over a 
large Tetf interval requires mixing to involve a given mass, ap- 
proximately the same independent of metallicity and Teff. 

4.3. Intermediate metallicity globular clusters 

For intermediate metallicity clusters, comparisons are made in 
Fig. [To] to patterns of anomalies for two stars in clusters[3 (M 13 
and NGC 1904) and in Fig.[T2]for observations of Fe as a func- 
tion of Teff in M 3, M 13, NGC 288 and NGC 2898. 

The expected patterns at six HB ages in a 0.55 Mq model 
with original [Fe/H] - -1.3 are compared, in the l eft panel o f 
Fig.[TOl with WF4-3085 from the cluster M13 ( BeM l2003h . 
According to Table 1 of that paper, the cluster iron abundance 
is [Fe/H] = -1.54. In the right panel, comparisons are also 
made to models with original [Fe/H] = -1.3 with star 469 of 
the cluster NGC1904 (M79) (Fabbian et al. 2005) which has 
[Fe/H] = -1.59 according to Kraft & Ivansi (^2003). Since ac- 
cording to Fig. [5] results are not too sensitive to metallicity (see 
also Fig.[T3ll this comparison is accurate enough given that mod- 
els with original [Fe/H] = -1.55 were not available. In WF4- 
3085 the abundances are well reproduced by models of 10 Myr 
or more. In 469, they are well reproduced by models of 5 Myr 
or more. Given error bars, all species heavier than Mg are well 
reproduced in both. One may argue about P in WF4-3085 which 
is claimed to be overabundant by a factor of 500 while the model 
gives a factor of 50 overabundance. The main difliculty is how- 
ever with He which is observed to be underabundant by a much 
larger factor than models predict. On the one hand this obser- 
vation is probably difficult since these are relatively cool stars 



* Variations in the abundances of Na and O, but not of Fe, as a func- 
tion of Tgg have recently been observed on the HB of M4 (Ma rino et al.l 
[2010). However this is for cooler stars than the Fe abundance varia- 
tions discussed here and in M4, the O and Na abundance variations are 
probably due to nuclear reactions. It would be of interest to verify the 
rotation velocities of the hotter stars in M4 given the correlations seen 
in Fig.\n\ 



where the He abundance is more difficult to determine. On the 
other hand. He is largey neutral in the atmosphere and this may 
enhance additional effets of diffusion there. These will be further 
discussed in section lSTI 

In the right panel of Fig. [12] the models with [Fe/H] = -0.7 
and [ Fe/H] =-13 are compared to observations of NGC 2808 
from Paceetal. (2006). The [Fe/H] = -1.14 determined by 
Can^etta et al. (2004) agrees within error bars with the [Fe /H] of 
the cooler HB stars (below 12000K) and is bracketed by tracks 
for models with [Fe/H] - 1.3 and -0.7. By zooming in on the 
figure, the origin (zero age HB) of each track may be clearly 
seen. The break between HB stars with the same Fe abundance 
as giants and those with an overabundance by a factor of ~ 20 
occurs slightly above 12000K which is ~ lOOOK hotter than 
in the three clusters M3, M13 and NGC288 in the left panel 
whose original [Fe/H] are bracketed by tra cks for [Fe/H] = -2.3 
and -1.3. In Fig. 4 of Ouiev v et al.l (l2009h . the Teft interval from 
10000 to 12000K is the interval over which anomalies are ex- 
pected in only a fraction of the stars. In M3, M13 and NGC288 
(left panel), all stars with T^ff > 1 lOOOK rotate slowly and have 
[Fe/H] within the range expected over most of the life time of 
a cluster with the original metallicity of the cluster The cooler 
stars have large enough rotation velocities for no anomalies to 
be expectecQ. Since we do not know the rotation velocities of the 
NGC 2808 stars between 11000 and 12000K, we do not know 
if the lOOOK difference in the Teff of the break in abundances 
is significant. One may then say that the run of Fe abundances 
in BHB stars shown in Fig. [H] from 1 1000 K to 17000 K is weU 
represented by calculations with a single value for the mixed 
mass of about 10"^ Mq. 

Finally the advantage of using clusters over field stars may 
be inferred from the right panel of Fig.[TT] There, the stars with 
Teff < 1 lOOOK give an indication of the original Fe distribution 
in the field. No information is however available on the origi- 
nal Fe abundance of individual stars with Teft > llOOOK. So 
the link wi th rotation is more d ifficult to establish. According 
to Fig. 4 of lOuievv et all (l2009l) . the 13,500 and 15000K stars 
with rotation velocities of 30^0 km s"' may or may not rotate 
fast enough for atomic diffusion to be affected, depending on 
whether meridional circulation penetrates or not in the surface 
convection zone. Since both stars might have a solar metallicity 
to start with, little can be inferred from their relatively large ro- 
tation velocities. The other field stars with large metallicity all 
have very small rotation velocities as was the case in clusters. 

4.4. Field sdB stars 

For sdB stars, neither the age nor the metallicity are known. 
Anomaly patterns are compared to observations for two stars 
with models of four different original metallicities. Comparisons 
are made with sdBs of Teff ~ 3 1000 K using data from 
O'Toole & Heber (2006) (Feige 48, see right panel of Fig.[5]l 
andTeff ~ 30000 K in Fig. [13] using data from Blanchet te et alJ 
(I2008h for PG0101-H039. We decided to compare to models of 
25 Myr, whose surface abundances are representative of a sig- 
nificant fraction of the evolutionary span. For Feige 48 (Fig.|5]l, 
14 of the 17 species which were both calculated (with origi- 
nal [Fe/H] between -0.7 and -1.3 or Zq between 0.001 and 
0.004) and observed agree reasonnably well. For PGOlOl -1-039 
(Fig.[T3ll. 15 species are both calculated and observed and 12 are 
in reasonable agreement with either of the two lower metallicity 



' For an interesting correlatio n between [F e/H], T^g and rotation in 
HB and sdB stars see Figure 8 of lCortes etaP ((2002). 
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Fig. 12 Left panel Surface concentration of Fe as a function of Teff in clusters (M3, M13 NGC288) from lBeh3l2003l with original 
Fe concentration between -2.3 and -1.3 dex the solar one. Models are shown for those two metallicities. Each color coded segment 
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Fig. 13 Surface abundances expected at 25 Myr for models with turb ulence anchored at AM/ Mg - 10 calculated with Zq = 
0.0001, 0.001, 0.004 and 0.02. Observations are for PGOlOl -1-039 from lBlanchette et al.1 (l2008l) . 



models. The agreement is however not perfect but this is perhaps 
not too surprising given the uncertainty of the age. 

Once neither age nor original metallicity are known, perhaps 
as useful information is obtained by comparing an ensemble of 
models spanning the metallicity, mass and age intervals of a 
large number of stars. Models of four metallicities are used with 
ages spanning the first 32 Myr of HE evolution and a mass in- 
terval leading to the observed rgft interval of 20000 to 37000 K. 
This Teff interval corresponds to that of a recent survey of the 
chemical composition of sdB stars whose preliminary results are 



published in 'Geie r et all (l2008h and lGeier et"an (1201 0^ with ad- 
ditional private communications to us from Dr Geier On Fig. [14] 
observations of X(Fe) are the grey (or red) circles with triangles 
being upper limits. The long dashed grey lines define the inter- 
val where observations are found. They are similarly defined for 
other elements from the observations and are shown in order to 
constrain the observed abundance interval in Fig.[T5]and[T6] The 
rotation velocities are shown on the lower panel of Fig. [14] They 
are all smaller than lOkm/s which is too small to interfere with 
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Fig. 14 Surface concentrations of Fe as a function of T^ff in sdB stars. Models were calculated with original concentrations of 
[Fe/H] =-2.3,-1.3, -0.4 and 0.0. Models span the T^g interval from 20 000 to 37 000 K and the mass interval from 0.47 to 0.52 Mq. 
The lines give the surface abundance variation for each calculated model. They are color coded as a function of the mass of the 
model which may be read by zooming on the figure in the electronic version; solid lines for the time interval from 10 to 32 Myr after 
ZAHB and dotted from to 10 Myr The dotted line starts at the original abundance of the model. For most models the time interval 
spanned is appro ximately 32 Myr. Note that most m odels starting with a solar metallicity have very short tracks. Observations 
(circles) are from iGeier et al.ll2008l iGeier et aDl2010l and Geier2010, private communication. Inverted triangles are upper limits. 
The stars in red in Fig.|2]are also in red in this figure. The bottom part of the figure gives the rotation velocities and the error bars on 



atomic diffusion at those T^fi according to lOuievv et al.l (l2009h . 
The uncertainty on T^-ff is also shown in that panel. 

Compari sons are made as a f unction of T ^ ff usin g observa- 
tions from IGeier et al.l (l2008l) and iGeier et all (1201 Ol) in Fig.[T5] 
for Ne, Mg, Al, Si, P, S, Ar, K, Ca,Ti, Cr and Fe and in Fig.[T6l 
for He and CN O. These are the s pecies for which observations 
are reported by iGeier et aL and which are also present in 

the OPAL opacity calculations. An indication of the limits of the 
domain of observed points is given by the long dashed gray lines 
which are defined for Fe in Fig. [14] For most atomic species, the 
lower broken gray line corresponds to as many upper limits as 
data points. We refer the reader to the original publication of the 
data. 

In all three figures, for each model and each chemical species 
the abundance is indicated by a curve originating at the original 
abundance, dotted for the first 10 Myr and solid until the end 
of the calculation, usually around 32 Myr The various colors 
distinguish masses (see footnote|5]l. The masses are indicated at 
the end of each colored curve on Fig.[T4]and may be read more 
easily by zooming in the electronic version. In doing the com- 
parison between observations and calculations, one should take 
into account that the calculations covered the first 32 Myr of 
the HB evolution. Since HB evolution probably lasts three times 



longer the colored lines should often be extended by a factor of 
2, or perhaps 3 when no saturation is involved. When saturation 
is important the lines can be much shorter however as is seen for 
many models which originally had the solar abundance (these 
are easily identified since their dotted line starts at the sohd hor- 
izontal gray line in Fig.fTSll. 

For Fe (Fig.[T4li. the calculated values cover the range of ob- 
served Fe abundances at all Teff's especially when the expected 
extension of the curves to the end of the HB evolution is taken 
into account. The concentration of Fe abundances between solar 
and 10 times below solar is to be expected given the expected 
original distribution of Fe abundances as seen in the right panel 
of Fig.[TT] A number of upper limits have been measured (the 
grey triangles) and can be accounted for by the early evolution 
of stars starting with very small original Fe abundances. Even 
those, however, end up spending most of their HB evolution 
within the main interval of the observations. The objects rep- 
resented by red circles are also compatible with the results pre- 
sented here even though they have a smaller gravity, and so are 
probably more evolved than the models (see Fig.|2]l. 

The comparison between calculated and observed abun- 
dances for elements from Ne to Fe is shown on Fig. [15] Iron is re- 
peated to allow evaluating the effect of representing observations 
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Fig. 15 Surface abundances for atomic species from Ne to Fe. Iron is also on Fig.[T4]where is found a detailed caption. The domain 
of observed values is bounded by long dashed broken gray lines which are also shown for Fe on Fig.[T4l Th e domain of abundances 
often extends below the lower curve since many measurements are upper limits (see iGeier et al. .2008 and lOeier et aDl2010.) . The 
solid horizontal grey line is the solar abundance. Note that individual panels span 7 orders of magnitude in abundances. 




Fig. 16 Surface concentration of He as a function of T^m in sdB st ars. See the caption of Fig.[T4]for the definition of the curves. 
Observations (circles) are from i Geier et al.l^OOSi iGeier et al? '2010 and Geier (2010, private communication). The stars in red in 
Fig.|2]are also in red in this figure. The lines for He were stopped between X(He) - 0.001 and 0.0001 because the evaluation of the 
surface He abundance becomes inaccurate for smaller values . 
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as lying between the two long dashed gray lines of each panel. 
However the lower bound frequently appears to correspond to 
the smallest abundance that can be detected since it generally 
correspo nds to uppe r limit s by the data po i nts wh ich may be 
found in iGeier etalj (l2008h and iGeier et al.l (|2010|) . The reader 
will find useful looking at the original data. In two cases. Si and 
S, observed abundances cover the same interval as the simula- 
tions' starting abundances (the dotted lines originate within or 
very close to the gray lines). In both cases they remain within 
those bounds. For Si, the small observed Teff dependence is ob- 
tained while for S, the observed slight increase around 35000 K 
is not reproduced. In four cases (K, Ca, Ti and Cr), the observed 
abundances are 4 to 5 orders of magnitude larger than the origi- 
nal abundances in the more metal poor original models. Taking 
into consideration the number of upper limits in the observa- 
tions, and that only the first third of HB evolution is covered by 
the lines, the agreement is generally good. One property which is 
not reproduced is the apparent increase in observed abundance 
above 33000 K for K and Cr For P and Ar the observed val- 
ues are compatible with observations except for Teft > 33000 K; 
however for P there are mainly upper limits. 

For C and N, (see Fig.fTSIl there seems to be strong disagree- 
ment for Teff > 33000 K while at lower Tetj, there is agreement 
for C while N is on average perhaps 3x more abundant than ex- 
pected. This could be related to some nuclear effects. On the 
other hand calculated O seems generally compatible with obser- 
vations. 

Our calculations of the surface He abundance compared to 
observations on Fig. [16] are constrained by a limitation of our 
code since in our calculation algorithm, X(He) is obtained by 
subtractioifl Consequently the lines for He were stopped be- 
tween X(He) = 0.001 and 0.0001 because the evaluation of 
the surface He abundance becomes inaccurate for smaller val- 
ues, when Y «. Z. One has the surprising result that observa- 
tions of X(He) are in agreement with calculations in the stars 
with Tgff > 29000 K but not cooler The agreement for stars with 
Teff ~ 30000 K is reinforced by the detailed comparisons shown 
on Figs.|5]and[T3] In the higher Teff stars, the observed abundance 
range is easily covered once one takes into account that calcu- 
lations cover only the first third of the HB life span. Whilst still 
helping at lower Teff, this is not sufficient to explain the absence 
of stars with an underabundance smaller than a factor of 0. 1 nor 
the presence of stars underabundant by a factor of 10 This re- 
inforces the problem with He abundance in globular cluster stars 
(see Fig.fTOll and is discussed in Sect. l5.ll 



5. Conclusion and General Discussion 

Stellar models have been evolved over the first 32Myrs of HB 
evolution for masses leading to the Teff interval of 8000 to 
37000 K and for metallicities ranging from Zo = 10""* to 0.02. 
A total of some 60 models were evolved. They all started from 
the evolution, from the zero age main-sequence to the He flash, 
of 0.8 to 1.0 Mq models of metalHcities from Zq = 10^"^ to 0.02 
described in iMichaud et al.l (12010 ) (see in particular their Sect. 
2). They corrrespond to the metallicities of the various globu- 
lar clusters whose horizontal branch has been studied. They also 
cover the metallicity (right panel of Fig.fTTTi and Teff intervals 
(right panel of Fig.|2]i of field sdB stars. 



^ It is the only element for which this is done. Eliminating this hmita- 
tion requires rewriting the central part of the evolutionary code involv- 
ing the transformation of H into He. 



5.1. The main results 

One remarkable observational property of both field sdB stars 
(studied here up to 37000 K) and cluster HB stars with T^ff > 
llOOOK is that their Fe abundance is relatively close to solar 
whatever the metallicity of the cluster and independent of the un- 
known original Fe abundance in the case of sdB stars. In the cal- 
culations described above for field stars, the original metallicity 
was varied from Zq = 0.0001 to Zq - 0.02. In globular clusters, 
the original abundances for the calculations are those of giants of 
the cluster. So long as HB stars have their outer 10^^ Mq mixed, 
it has been shown to follow naturally from the stellar evolution 
with atomic diffusion that the final Fe abundance ended within 
a factor of ten interval nearly independent of Teff (see Sec. |4] and 
in particular Figures [TT] [12] and [14). Overabundances of Fe by 
factors of up to 100 are implied. Reducing the mixed mass by a 
factor of ten or more leads to unacceptably large Fe abundances 
(see Sec. 13. 21 and Fig.|6]l. The well observed Fe abundance so 
constrains the mixed zone in stellar evolution modelfl Is this 
mixing compatible with observation of other atomic species? 

Detailed comparisons were carried out for a number of sdBs 
and cluster stars in which the abundance of many elements were 
determined. In all cases the comparison between observations 
and model calculations led to acceptable results (see Sec. |4] and 
Figs.l5ll9] [To] and [T3Tl. Approximately 12 of 15 elements are in 
reasonable agreement in each case. This is obtained without any 
adjustment since the only adjustable parameter was fixed by the 
Fe abundance. This is true in globular clusters whose original 
[Fe/H] varies from -2.3 to -0.7 as well as in field sdB stars. For 
these, a comparison to a large sample of field stars covering a 
large age and original metallicity interval was also carried out 
with general agreement (see Figs. [15] and [16). 

The calculated surface abundance of He agrees with that ob- 
served in the sdB stars shown in Figs. [5] and [13] and also for 
Teff > 25000 K in Fig.[T6]but there appears to be a large discrep- 
ancy in the cooler stars and in particular those HB stars shown in 
Fig.[TO]and the left panel of Fig. [5] In the latter cases, Teff is small 
enough that He is not ionized in the atmosphere. The neutral He 
atom ic diffusion co efficient is ~ xlOO larger than the ionized 
one (Michau d et alj|197 8). Additional element separation could 
be occuring in the stellar atmosphere for He since its settling ve- 
locity there would be enhanced by its being neutr al. There is ob- 
servational evidence for such stratification of Fe (iKhalack et al.l 
2007, 2010; LeBlanc et al. 2009, 2010). Other species, for which 
grad in the atmosphere would happen to be very large, could also 
be affected by separation in the atmosphere if the mixing there 
is not very strong. The evolutionary models described in this pa- 
per may then not tell the whole story. In real stars, there could 
be additional abundance variations between the surface and the 
bottom of the mixed mass of some 10"^ Mq. 

Since in particular gravity varies by two orders of magnitude 
over that Tetj interval on the HB, the similarity of observed abun- 
dances is quite surprising. It is clearly linked to the saturation of 
radiative accelerations. 



' If, instead of using Eq. ((3) one used Eq. {2) so that turbulence were 
fixed at a given To, the mixed mass would vary by more than 3 orders 
of magnitude between stars of Tgff = 12000 and 30000 K. We have 
verified that starting with Zq = 0.001, a 12000 K HB star would develop 
solar surface Fe abundance while a 30000 K star would develop a Fe 
abundance 50 times solar. This would lead to completely unacceptable 
changes in relative surface abundances when taking sdBs and HB stars 
of clusters into account (see Figs. 1 121 and [T4l for observed surface Fe 
abundances and also Fig.|6]for the effect of a change of the mixed mass 
by 2.5 orders of magnitude). 
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The agreement between the observations over the whole HB 
and the expected abundance anomalies from the evolution cal- 
culations leaves little doubt that radiative accelerations are the 
main cause of abundance anomalies on the HB. 

5.2. Turbulence or mass loss 

In this paper it was found that if the outer 10"^ Mq was mixec(3 
by turbulence most observed abundance anomalies followed. 
However as soon as an adjustable parameter has an effect on the 
results, one must remain cautious and, in spite of the observa- 
tional support described above, one must question if the physical 
process assumed to do the mixing is the only one possible or if 
there exist alternatives. In the only stellar evolution calculations 
we know of where mixing (or homogeneity) was n ot assumed 
for the outer region, it has been shown by IVick et al. (2010) 
that abundance anomalies observed on AmFm stars could be re- 
produced as well by a model that assumes mass loss as by the 
model where turbulence is in compet ition with atomic diffusion 
([Richer etal .1120001: iTalon et al.1 12006). They found that the ve- 
locity corresponding to the mass loss rate had to equal the atomic 
diffusion velocity at the stellar mass fraction where the mixed 
zone ends in the turbulence model. This occurs here at 10"^ Mq. 
Presumably the same would be true in HB stars. Instead of tur- 
bulence, mass loss could play a role to reduce anomalies. One 
may evaluate the required mass loss rate by equating the corre- 
sponding velocity to the He settling velocity at IQ^^ Mq using 
the expression 



— = -47rr pvdrift(He). 
at 



(4) 



It was found to equal 3-5 10"''*MQ/yr at log AM/M, = -7, in 
a model with Teft ~ 30000 K. A ten times smaller mass loss 
rate has been suggested before to explain in p articular observa- 
tions of Si in sdB stars (iMichaud et al.l[T985l) . The calculations 
included only the outer 10"^ Mq and the separation occured right 
in the atmosphere. Both the mass loss rates and the timescales of 
the calculations were smaller A complete stellar evolution cal- 
culation including mass loss instead of turbulence as the mixing 
process seems justified for HB stars, along the lines done by 
IVicketal.1 (120101) for main-sequence stars. 



5.3. fi gradient inversions 

It was shown in Sect. l3.2l that when the mixed mass is smaller 
than 10^^ Mq an inversion of the mean molecular weight often 
oc curs (see F ig.[6]l. It has been suggested in a similar context 
by iTheado et al. (2009) that such a fio inversion is unstable and 
should lead to mixing. It is tempting to conclude that this /iq 
inversion is the cause of the mixing implied by observations. 
However their analysis does not include the effect of g^d on the 
instability. Since the radiative accelerations are the cause of the 
fio inversion, their role in the analysis of the hydrodynamical in- 
stability could be important. The metals which are supported by 
gradS do not Contribute to the increase of the weight of the ma- 
terial. So one may wish to consider an effective fio for which 
supported metals contribute negatively. Intuitively, one may as- 
sume that if the function 



/^Oeff 



Y,N{AdAimp[g-gM)] 
Y,N{Ai)mpg 



(5) 



does not decrease inwards where yuo decreases inwards because 
of giadS (which cause some metal abundances to increase out- 
wards) then gradS may maintain stability. This function yUoes is 
shown (inverted) by the dotted curves (see Fig.|6]l. Contrary to 
1 I^Q, 1 / /loeff increases outwards where the Fe and Ni abundances 
increase outward. It seems to us that one should not conclude 
that a inversion caused by gndS on metals leads immediately 
to an instability; further work is needed to investigate the nature 
of the ju inversion instability in the presence of radiative acceler- 
ations. 

5.4. Asterosismology 



The first mesh point in the models of this paper was typically at 
10-'^ Me. 



In their ground breaking work, ICharpinet et al.l (Il996b have es- 
tablished that the Fe abundance expected to follow from atomic 
diffusion in the envelope of sdB stars should le ad to pulsations. 
This was confirmed observationnally ( Kilkenny et al. 199 7|). See 
also th e reviews of ICharpinet et al.l (12001 1) and .Fontaine et al.l 
(l2008h . Their approach has the advantage to involve no ad- 
justable parameter since they assume equilibrium between ^rad 
and gravity in the driving region for the pulsations. This im- 
plies that mass loss and turbulence have a negligible effect. It 
is however not established if in sdBs with small original metal- 
licity there is enough Fe to fill sufficiently the region where 
the driving occurs. As noted in the conclusion to the tenth 
sdB f or which a detailed asterosismic analysis had yielded the 
mass dRandall et alj2009l) a more accurate fit to the asterosismic 
data requires evolutionary models of sdB stars. The appr oach to 
equi librium in the driving region was studied by Fontain e"et al.l 
(2006) starting from a solar Fe abundance. They showed using 
static models that the Fe abundance became able to drive pulsa- 
tions within 10^ yrs. These works are based on a detailed analy- 
sis of the structure of the outer regions of these stars. 

Our results cleary show that from whatever Fe abundance 
sdB stars start with, a sdB star model with the outer 10"^ Mq 
mixed ends up with a solar Fe abundance in a fraction of the HB 
lifetime of 100 Myr T his could be used as a starting abundance 
for the calculations o flFontai ne et al.'(2006') mentioned above. In 
fact, the curve labeled -9 in Fig.|6]shows that after 32 Myr, start- 
ing from a metallicity 20 times below solar, the complete evolu- 
tionary models calculated here with OPAL data have accumu- 
la ted approximately the same Fe abundance as shown in Fig. 1 
of Fontaine et al. (2006) at the same age for a model starting with 
a solar abundance. According to Fig. [14] models starting with an 
Fe abundance 200 times below solar end up wi th 3 times less 
Fe afte r 32 Myr but according to Fig. 1 and 3 of lFontaine et al.l 
(120061) this is still sufficient to excite the spectrum. 

More recently, using an analysis of the F e abundance re- 
quired to cause observed pulsations in sdBs, ICharpinet et al.l 
(l2009l) conclude that if Fe is the sole atomic species responsi- 
ble for the pulsation of sdB stars, the Fe abundance needs to be 
log%g/A^H ~ -4.09 and -3.75 at Teff = 29580 and 35050 K re- 
spectively in the driving region. This is larger than the observed 
Fe abundance and suggests that some separation needs to occur 
between the driving region and the surface. However from Fig .[8] 
Ni is easily a factor of 3 more overabundant than Fe and it may 
contribute significantly to driving the pulsations as suggested by 
iJefferv & Said (l2007h and just as it contributes s ignificantly to 
iron c onvection zones in Pop. I stars (see Fig. 4 of iRichard et al.l 
120011) . 

The models described in this paper clearly need to be tested 
by calculating their pulsation spectra. When models with mass 
loss become available (see Sect. l5.2l i. asterosismology could 
help decide which is the main mechanism competing with 
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atomic diffusion driven by brad's. This is however outside the 
scope of the present paper. 
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Appendix A: Black and white version of FigureS] 
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Fig. A. 1 Black and white version of Fig.|4] Gray coded concentrations in two HB stars of the same metallicity after 25 Myr on the 
HB. Left panel with a T^sf of 14 000 K (0.59 Mq) and right panel of 30 100 K (0.51 Mq). The scale of the radius is Hnear but the 
logarithmic value of the mass coordinate above a number of points, log AM/M*, is shown on the left of the horizontal black line. The 
concentration scale is given in the right insert. Small circles near the top of the left panel mark the extent of the surface convection 
zone while similar circles near the center of both models mark the central convection zone. The small inset in between the two 
panels shows the high Teff star on the radius scale of the low T^ff star. For -7 < log AM/M, < -4 the concentration is quite different 
for many species. It is surprisingly so for C and O for log AM/M, > -2. 
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Appendix B: Properties of evolutionary models 

Some of the properties of all evolutionary models calculated, 
grouped by metallicity. All models in Figs. IB. 1 I to IB .41 were cal- 
culated with a turbulent diffusion coefficient given by equations 
O and (O with AMq = 10"^^ Mq. See Sect.i] 
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Table B.l. Models illustrated in this article, with their ZAHB characteristics 
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log* 
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0.0001 
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10036 


3.601 




13.194436 


209248 


0.0001 


0.6900 


45.115 


9050 


3.404 




13.194433 


209251 


0.0001 


0.7300 


50.354 


7 480 


3.050 




13.194448 


209252 


0.001 


0.497455 


19.181 


36542 


6.058 


dM-7.5D10K-4 


12.204222 


209382 


0.001 


0.49746 


19.179 


36 078 


6.036 


" 


12.204233 


209322 


0.001 


0.4975 


19.184 


35057 


5.986 


" 


12.204219 


209314 


0.001 


0.4980 


19.250 


32405 


5.848 




12.204192 


209313 


0.001 


0.4982 


19.276 


31898 


5.820 




12.204193 


209312 


0.001 


0.4985 


19.309 


31288 


5.786 


" 


12.204252 


207669 


0.001 


0.4990 


19.369 


30494 


5.741 


" 


12.204227 


205714 


0.001 


0.5000 


19.482 


29 310 


5.670 


dM-lO.ODlOK-4 


12.204068 


209839 


0.001 


0.5000 


19.482 


29318 


5.671 


dM-9.0D10K-4 


12.204068 


209840 


0.001 


0.5000 


19.482 


29322 


5.671 


dM-8.0D10K-4 


12.204068 


209843 


0.001 


0.5000 


19.478 


29324 


5.671 


dM-7.5D10K-4 


12.204182 


207709 


0.001 


0.5025 


19.753 


27 377 


5.548 




12.203995 


211226 


0.001 


0.5050 


19.984 


26007 


5.456 




12.203976 


208085 


0.001 


0.5100 


20.418 


24005 


5.312 




12.203820 


207668 


0.001 


0.5200 


21.238 


21 284 


5.094 




12.203564 


207704 


0.001 


0.5300 


22.135 


19 343 


4.918 




12.203353 


207706 


0.001 


0.5500 


24.927 


16403 


4.596 




12.203071 


207707 


0.001 


0.5600 


27.324 


15043 


4.414 




12.202989 


207467 


0.001 


0.5700 


30.435 


13627 


4.203 




12.202952 


205899 


0.001 


0.5800 


33.931 


12169 


3.967 




12.202947 


206117 


0.001 


0.5900 


37.393 


10734 


3.714 




12.202944 


207665 


0.001 


0.6000 


40.585 


9 387 


3.453 




12.202956 


208081 


0.001 


0.6100 


43.386 


8 166 


3.189 




12.202966 


208083 


0.004 


0.48090 


16.254 


35 068 


6.044 


dM-7.5D10K-4 


11.072471 


209535 


0.004 


0.48095 


16.244 


33 257 


5.952 




1 1 .07267 1 


206918 


0.004 


0.4810 


16.258 


32667 


5.921 




11.072670 


206852 


0.004 


0.4820 


16.368 


29 282 


5.728 




11.072520 


206850 


0.004 


0.4835 


16.516 


27 311 


5.605 




11.072435 


211224 


0.004 


0.4850 


16.647 


26021 


5.519 




11.072454 


206849 


0.004 


0.4875 


16.854 


24447 


5.407 




1 1 .072200 


211225 


0.004 


0.5000 


17.805 


19 923 


5.039 




11.071824 


206917 


0.004 


0.5200 


20.075 


15 842 


4.606 




11.071256 


206832 


0.004 


0.5400 


28.002 


11 683 


3.949 




11.071092 


206914 


0.004 


0.5600 


36.861 


7 265 


3.020 




11.071109 


206960 


0.004 


0.5700 


39.702 


5 793 


2.602 


„ 


11.071069 


206962 


0.02 


0.472895 


14.060 


32783 


5.982 


ciM-7.5D10K-4 


11.712004 


209540 


0.02 


0.472897 


14.061 


32547 


5.970 




11.712016 


209539 


0.02 


0.47290 


14.064 


32347 


5.959 




11.712347 


203666 


0.02 


0.47293 


14.061 


31429 


5.909 




11.711976 


204568 


0.02 


0.47295 


14.065 


31 108 


5.891 




11.711949 


200008 


0.02 


0.4730 


14.070 


30552 


5.860 




11.711964 


199964 


0.02 


0.4750 


14.253 


25 862 


5.566 




11.711779 


199984 


0.02 


0.4900 


15.363 


18 570 


4.972 




11.711468 


204456 


0.02 


0.4950 


15.781 


17 268 


4.838 




11.710998 


199777 


0.02 


0.5260 


30.576 


8708 


3.388 




11.710788 


204509 



*Zo is the pre-main-sequence starting metallicity of the models, before applying or-corrections. Note that 
no such correction is applied to Zq = 0.02 models. 

''From zero age main-sequence. 

'^Identification number which is used to identify models on some of the figures. 
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Fig. B. 1 Models with Zq = 0.0001 . The curves are identified on the figure by the mass of the HB model. The color code is the same 
as used for the left panel of Fig.|2]and is defined in footnote ||5|. Zero age HB is defined as the moment when X('^C) reaches 0.04. 
On the figures the curves are not all very precisely aligned to that value. Xo(surf Fe) identifies the Fe mass fraction in the original 
main-sequence model. Similarly for the other species. 
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Fig. B.2 Models with Zq = 0.001. See caption ofFig.lBJl 
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Fig. B.3 Models with Zq = 0.004. See caption of Fig.lBH 
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. B.4 Models with Zq = 0.02. See caption ofFig.lBll 



